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Highly relativistic jets are most probably driven by strong magnetic fields and 
are launched from the accretion disk surrounding a central black hole. 

In this paper we review some of our recent results considering the two-dimensional 
magnetic field structure and the dynamics of collimating relativistic jets. Apply¬ 
ing the stationarity assumption enables us to calculate global solutions of the 
relevant MHD equations. The structure of the jet magnetosphere follows from 
the so called Grad-Shafranov equation for the force-balance between axisymmet- 
ric magnetic surfaces. The plasma dynamics along the magnetic field lines is 
given by the solution of the MHD wind equation. In the force-free assumption, 
which is appropriate for relativistic jets, both equations de-couple. 

We discuss solutions of the Grad-Shafranov equation obtained in a general rela¬ 
tivistic context applying the 3-1-1 formalism for Kerr geometry. These solutions 
extend from the inner light surface around the Kerr black hole to the asymptotic 
regime of a cylindrically collimated jet with a finite radius. 

In a further step, we include differential rotation of the foot points of the field 
lines. In this case, the shape of the light surface is not known a priori and must 
be calculated in an iterative way. The solutions show that differentially rotating 
jets are collimated to smaller radii compared to jets with rigid rotation. 
Gonsidering the general relativistic MHD wind equation, we investigate the dy¬ 
namics of the collimating jet, in particular the effects of Kerr metric on the 
acceleration. Temperature and density follow a power law versus radius. The jet 
magnetic (velocity) field is dominated by the toroidal (poloidal) component. 
Having at hand a relativistic MHD jet solution, we calculate the thermal optically 
thin X-ray spectrum for the innermost hot part of the jet. Doppler shift and 
boosting is taken into account. For microquasars we obtain a jet X-ray luminosity 
« 10^^ergs“^. Iron emission lines are clearly visible. 


1 Relativistic astrophysical jets 

Relativistic jets originating in the close environment of a rotating black hole 
are observationally indicated for two categories of sonrces concerning mass 
and energy ontpnt. In active galactic nuclei (AGN), the jets are launched 
in the magnetized environment around a rotating, super massive black hole 
(Blandford & Payne 1982, Sanders et ah 1989). The discovery of Galactic 
microquasars has proven superluminal motion on a much smaller energy scale 


(Mirabel & Rodriguez 1994), the central engine, however, is believed to work 
similarly to the AGN. 

From the observations we know that astrophysical jets are collimated to al¬ 
most cylindrical shape. Relativistic jets are generally detected in non-thermal 
radio emission, clearly indicating the magnetic character of jet formation - 
magnetic acceleration and collimation. There is clear evidence that jet for¬ 
mation is always connected to the existence of an accretion disk. As the 
central object in AGN and microquasars is a black hole, the surrounding ac¬ 
cretion disk is the only possible location for a magnetic held held generation. 
Gonsidering the before mentioned observational constraints, it is clear that a 
theoretical, quantitative analysis of the jet structure in these sources must take 
into account both magnetohydrodynamic (MHD) ehects and (general) relativ¬ 
ity. Relativistic MHD implies that electric helds, which are unimportant in 
Newtonian MHD, cannot be neglected. 

Due to the complexity of the MHD equations, all stationary, relativistic 
models of magnetic jets so far, rely on simplifying assumptions such as self¬ 
similarity (Gontopoulos 1994), some other prescription of the held structure 
(Li 1993, Beskin 1997) or the restriction to asymptotic regimes (e.g. Appl 
& Gamenzind 1993; Nitta 1997; Beskin & Malyshkin 2000). As relativistic 
jets must be highly magnetized, the force-free limit of MHD may be applied. 
This allows for a truly two-dimensional (2D) calculation of the held structure 
(Gamenzind 1987; Fendt et al. 1995; Fendt 1997; Fendt & Memola 2001), in 
contrast to the self-similar treatment. 

The stationarity assumption itself is essential in order to obtain a global 
solution of the equations. Yet, time-dependent MHD simulations of relativistic 
jet formation are limited in spatial scale and to a short time evolution (Koide 
et al. 2000). On the other hand, we do not really know whether a stationary 
MHD solution obtained is actually stable. 

2 Stationary relativistic MHD 

The axisymmetric, stationary electromagnetic equations in the context of Kerr 
metric were hrst derived by Znajek (1977). The more comprehensible for¬ 
malism of a “3-1-1 split” has been introduced by Thorne et al. (see Thorne 
et al. 1986). Gamenzind (1986, 1987) formulated a fully relativistic station¬ 
ary description of axisymmetric MHD hows presenting solutions of the Grad- 
Shafranov and the wind equation. An essential step to obtain the global mag¬ 
netic held structure is to take into account properly the regularity condition 



along the light surface. This matching problem of relativistic magnetospheres 
has been solved by Fendt et ah (1995). The general relativistic version of 
the Grad-Shafranov equation including inertial terms and entropy has been 
presented by Beskin & Pariev (1993). Solutions of the wind equation in Kerr 
geometry considering the stationary plasma motion along the magnetic held 
were obtained by Takahashi et ah (1990), however, mainly discussing the ac¬ 
cretion how onto the black hole. 

2.1 Space-time around rotating black holes 

In the 3-1-1 split the space-time around rotating black holes (mass M, angu¬ 
lar momentum per unit mass a = J/Mc) can be described using the Boyer- 
Lindquist line element[] 

ds^ = a^c^dt^ — {d(j) — oodt^ — (p^/A) dr'^ — dO'^ . (1) 

t denotes a global time in which the system is stationary, (j) is the angle around 
the symmetry axis, and r, 6 are similar to the hat space spherical coordinates. 
The electromagnetic held B,E, the current density j, and the electric charge 
density pc are measured by the ZAMO according to the locally hat Minkowski 
space. These local experiments are put together by a global observer at a 
certain global time using the lapse and shift function to transform from local 
to global frames. In spite of this transformation. Maxwell’s equations in the 
3-1-1 split look similar to those in Minkowski space. 

2.2 The cross—field force—balance 

The axisymmetric held structure follows from the force-balance across mag¬ 
netic hux surfaces. The magnetic hux function \k(r, 0) = ^ J Bp-dA measures 
the magnetic hux through a loop of the Killing vector rh = Similarly, 

the total poloidal current is dehned by integrating the poloidal current density 
over the same loop I = — f ajp ■ dA = —^auBp. The indices P and T de¬ 
note the poloidal and toroidal components of a vector. In the force-free limit, 
PcE + X B = 0, we have jp || Bp and, thus, I = /(T). For a degenerated 

^ The parameters of the metric tensor are defined as usual, cos^ 6, A = 

- 2GMr/&‘ + = (r^ -p - a^A sin^ 6», w = (E/p)sin6l, w = 

2aGMr/cY?, a = p\fK/T,. Here, Lo is the angular velocity of an observer moving with 
zero angular momentum (ZAMO), uj = (dd/dt)zAMO and a the lapse function, describing 
the lapse of proper time r in the ZAMO system to the global time t, a = {dT/dt)zAMO- 



magnetosphere, \B‘^ — E‘^\ » \E ■ B\ ~0 another conserved quantity can be 
derived. The derivative of the time component of the vector potential dehnes 
the iso-rotation parameter hip = = —c{dAQ/d^), which is sometimes 

interpreted as angular velocity of the magnetic held lines. 

The force-balance across the held can be derived from the toroidal com¬ 
ponent of the stationary Ampere’s law and is commonly known as Grad- 
Shafmnov equation (GSE) for the hux function \h(r, 6), 

(hv ■ = (h I v^|2 _ 2 ^2) 

\ uj^ J aG aoo 

(see Okamoto 1992). Here, D = 1 — (cD/cDl)^ The ’ denotes derivatives d/dT, 
and a;L(T) the position of the two light surfaces Q, oj^ = (±a c/(Hp — ■ 

2.3 Jet acceleration — the force—balance along the field 

The stationary, polytropic, general relativistic MHD how along axisymmetric 
hux surfaces T can be described by the wind equatior^ 

fEykok2-2k2M^-hM^ ,, 

for the poloidal velocity Up = yUp/c (Camenzind 1986, Takahashi et ah 1990). 
The Alfven Mach number M is dehned via = Afinu^/B"^^ with the proper 
particle density n, the specihc enthalpy p, and the poloidal magnetic held 
Bp = B-p/{gQQ) -t-^fos^p), rescaled for mathematical convenience. The specihc 
total energy density E{'^) and total angular momentum density T(\k) are con¬ 
served along the hux surfaces. For a polytropic gas law (T = n/m), the wind 
equation can be converted into a polynomial equation ~ ^ 

(Camenzind 1986), where the coefficients Ai depend on the radius along T, 
the held strength and the parameters Qp, E, L,a^. The magnetization cr*(T) 
measures the Poynting hux in terms of mass hux. 

3 The 2D structure of relativistic jet magnetospheres 

We now discuss the GSE solution for a force-free, collimating jet in the general 
relativistic context (see Fendt 1997), assuming Gp’(T) =const. As asymptotic 

^ The ‘-I-’ (‘—’) sign holds for the outer (inner) light surface with ftp > ^ (resp. ftp < 

^ with the following abbreviations, fco = + 2 go 3 flF + ffoo: = 1 — ^viL/E), and 

ki = — (gas + 2gaa{L/E) + gooG/E)"^) / (^Qg — googas) 






0 12 3 



Figure 1: Axisymmetric magnetic field structure of relativistic jets '^{R,Z). 
General relativistic jet magnetosphere around a rotating black hole with a = 
0.8, Gp =const. {left). Special relativistic magnetosphere with differential 
rotation Gf('1/) {right). 


boundary condition we apply the analytical, special relativistic ID GSE solu¬ 
tion of Appl & Gamenzind (1993). Since / = /(d'), this asymptotic solution 
provides the GSE source term also for the collimation region. The disk mag¬ 
netic held distribution is parameterized as T ~ {x/b)'^/{l + {x/h)'^) + 
Here, the magnetic hux anchored in the black hole is = 0.5 with m = 3 
and a core radius b. The asymptotic jet radius can be parameterized in terms 
of the light cylinder radius Ri or the gravitational radius M. For a = 0.8, 
(Df/Dh) = 0.4, the jet radius is 3i?L corresponding to 30 M. 

The main features of the calculated jet magnetosphere are the following 
(Fig. 1, left). The held lines originate near the inner light surface close to 
the rotating black hole and collimate to an asymptotic jet of hnite radius of 
several (asymptotic) light cylinder radii. The solution is dehned on a global 
scale, satisfying the regularity condition along the light surfaces. We hnd a 
rapid held collimation within 20 M distance from the source. The near-disk 
solution has three diherent regimes. Here, the magnetic hux is either outgoing 
towards the asymptotic jet or in-going towards the black hole. But there exist 
also hux surfaces near the axis which not connected to the disk. Thus, if we 
imagine a mass how associated with the held lines, we expect three diherent 
how regimes - accretion, outhow, and a region empty of a plasma how. 




Extending our work on relativistic magnetospheres, we have investigated 
jets with non-constant (see Fendt & Memola 2001). Here, general 

relativity is not taken into account. Since astrophysical jets seem to be always 
connected to an accretion disk, differential rotation of the held lines should 
be a natural ingredient for any jet model. However, as a difficulty, then the 
position and shape of the light surface singularity is not known a priori, but 
have to be calculated iteratively with the held distribution. Similar to /(fh), 
also the rotation law HF('h) can be taken from the asymptotic GSE solution. 

The calculated magnetic held structure is shown in Fig. 1 {right). The 
diherentially rotating magnetosphere is collimated into a narrower jet. This, 
however, can be balanced by an increase of electric current. Comparison to 
high-resolution radio observations of the M87 jet (Junor et ah 1999) shows 
good agreement in the collimation distance. For M87 we derive a light cylinder 
(jet) radius of 50 (120) Schwarzschild radii. 

The held structure is governed by /(T) and r2F('h)- In combination with the 
disk magnetic hux distribution this allows to determine the magnetic angular 
momentum loss from the disk in the jet and the toroidal magnetic held along 
the disk. The angular momentum hux per unit time per unit radius is dJ/dx = 
—xBJ{x) along the disk. In our solutions, most of the magnetic angular 
momentum is lost in the outer part of the disk (Fendt & Memola 2001). 

3.1 The MHD wind solution in Kerr metric 

Here, we discuss solutions of the stationary, general relativistic, MHD wind 
equation (3) along collimating magnetic hux surfaces. The solutions are calcu¬ 
lated on a global scale up to several 1000 gravitational radii. Diherent magnetic 
held geometries were investigated, parameterized by the shape of the magnetic 
hux surface and the hux distribution (Fendt & Greiner 2001, FGOl). 

Applying microquasar parameters, we obtain the following results (Fig. 2). 
The solution passes all three critical points (slow, Alfven and fast magnetosonic 
point). Substantial acceleration is achieved also beyond the Alfven point due 
to the magnetic nozzle ehect (Gamenzind 1986, Li 1993). The temperature in 
the inner jet is up to lO^^K. The jet injection radius is at 8.3 M. Gomparison 
to solutions with smaller injection radius (3.3 M) allows to clarify the role of 
general relativity for the jet acceleration. In this case, the asymptotic velocity 
is much higher due to the faster rotation Df at smaller radii. With Df = 0.14 
the asymptotic velocity increases to tip = 8.48 (FGOl). 

In the limit of Minkowski metric the jet magnetization can be low, although 
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Figure 2: MHD wind solution in Kerr metric (a = 0.8). Shown is the normal¬ 
ized poloidal velocity Wp in the global frame {left) and the density {thick line) 
and temperature versus radius x = r/M {right). The intersections of the phys¬ 
ical wind branch with the secondary solution indicate the magnetosonic points. 
Dynamical parameters are (see FGOl); n* = 980, Dp = 0.035, injection sound 
speed Cs* = 0.05, Alfven radius xa = 22.9, asymptotic velocity Up^o = 2.53, 
total energy E/rripC^ = 2.79, held distribution ~ and Bp ~ a; 


the asymptotic speed is the same as for Schwarzschild metric (FGOl). The jet 
how is not ahected by gravity and, thus, needs less magnetic energy to gain 
the same asymptotic speed by magnetic acceleration. Solutions for diherent 
Kerr parameter a show that the jet is faster for smaller a (see FGOl). The 
reason for this seems to be the fact that the ehective potential of a black hole 
weakens for increasing values of a. 

These results might not be surprising as just demonstrating the fact of 
a magnetically driven jet. MHD theory tells us that magnetic acceleration 
takes place mainly around the Alfven point. If that is far out, the inhuence of 
general relativity must be marginal. 

3.2 X-rays as tracer for jet motion in microquasars? 

The MHD wind solutions discussed above provide the how dynamics along 
a prescribed poloidal magnetic held line with maximum temperatures up to 
10^^ K in the innermost part of the jet. Here, we calculate the thermal optically 
thin X-ray spectrum of such a jet. Prescribing the jet mass how rate Mj 
together with the shape of the held line, the wind solution gives a unique 









Figure 3: X-ray spectrum of a relativistic MHD wind solution for a micro¬ 
quasar. Shown is the luminosity of two volume elements along the jet with 
different temperature, T = 10^K {left) and T = 10®K {middle), and veloc¬ 
ity for different inclination, i = 40°, 20°, 0°,—20°,—40° {dashed cnrves from 
top to bottom), and in the rest-frame {thick curve). The right panel shows 
the overall spectrnm of the red (blue) shifted side of the jet, i = — 20 ° {top) 
and i = — 20 ° {bottom), the combined spectrum (thick line middle) and the 
rest-frame spectrum (thin line middle). 


set of parameters of the flow. For each volnme element of the expanding 
jet flow with decreasing temperature and density and increasing velocity we 
compnte the continunm spectrum and the emission lines of an optically thin 
plasma (see Memola et ah 2002). The spectra of the single volumes are then 
put together to an overall spectrum taking into acconnt the Doppler factor 
D = ( 7(1 — /3cos9))~^, for shifting {Eo = D Ef.) and boosting {Lo{Eo) = 

Le{Ee)) the emitted energy Eg and luminosity Le- These factors vary dne 
to the different velocity (3 of the matter in these volnmes and the different 
inclination 9 of the velocity vectors to the line of sight. 

Figure 3 shows the spectra of two example volume elements of different 
temperatnre (and velocity) for different jet inclinations. The combined total 
spectrum of all volume elements within the collimating jet cone is shown in 
Fig. 3 {right). We hnd a X-ray luminosity ~ 10^^ (Mj/lO”® M 0 yr“^)ergs“^ for 
a microquasar. The 6.6 and 6.9 keV emission lines can be identihed as Ka 
iron lines, while the one at 8.2 keV could be the Kft. We emphasize that our 
approach is not a fit to certain observed spectra. In contrary, for the hrst time, 
for a jet flow with characteristics dehned by the solntion of the MHD wind 
equation, we derive its X-ray spectrum. 
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